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Abstract 

In this paper we construct a class of random matrix ensembles la- 
belled by a real parameter a G (0,1), whose eigenvalue density near 
zero behaves like 1x1°^. The eigenvalue spacing near zero scales like 
and thus these ensembles are representatives of a continu- 
ous series of new universality classes. We study these ensembles both 
in the bulk and on the scale of eigenvalue spacing. In the former case 
we obtain formulas for the eigenvalue density, while in the latter case 
we obtain approximate expressions for the scaling functions in the mi- 
croscopic limit using a very simple approximate method based on the 
location of zeroes of orthogonal polynomials. 



1 Introduction 

Random matrix ensembles arise in an overwhelming number of diverse appli- 
cations. Their utility stems from the fact that even though one can construct 
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a multitude of different random matrix models, when one studies properties 
on the scale of eigenvalue spacing the results become universal and indepen- 
dent of the detailed structure of the random matrix measure. Hence one 
can extract predictions even without knowing the details of the appropriate 
microscopic model of the phenomenon being studied. 

The type of universality class depends in general on the symmetry prop- 
erties of the random matrix model and, of course, on the scaling proper- 
ties of eigenvalue spacing. Thus one has different bulk and edge universal- 
ity regimes, as well as a discrete series of multicritical universality regimes. 
Consequently only a discrete series of (rational) scaling exponents could be 
obtained. 

In this paper we would like to show that there exists a real continuous 
range of universality regimes with generic real scaling exponents. We will 
perform an explicit construction for ensembles with eigenvalue spacing be- 
having like l/iV^/(i+°) with a E (0, 1). 

A similar continuous series of universality classes appeared naturally in 
for random matrix models with noncompact eigenvalue support and 
power-law tails in the eigenvalue density. There the universal behaviour ap- 
peared for large eigenvalues. In general, however, the appropriate measures 
were difficult to obtain in an explicit way. This was the main motivation for 
this investigation. 

In this paper we would like to perform a construction of a more conven- 
tional class of models with compact support and where the universal regime 
is located close to A = 0. 

We strongly suspect that a mapping of this behaviour to infinity would 
correspond to the universal scaling for the Levy ensembles of [|1|. Another 
possible application of similar random matrix models (in their chiral variant) 
might be to model the behaviour of low lying eigenvalues of the QCD Dirac 
operator at the chiral phase transition. In this context the exponent a in 
p(A) ~ I A I" is the inverse of the S critical exponent of the chiral phase 
transition [|, |, |. 

The plan of this paper is as follows. In section 2 we briefiy recall the main 
features of the old and new universality regimes in random matrix theory. In 
section 3 we construct the random matrix measures of the new multicritical 
ensembles and compute the eigenvalue density in the bulk of the spectrum. In 
section 4 we move to a discussion of the microscopic (presumably universal) 
scaling behaviour and propose a simple approximate method based on the 
distribution of zeroes of orthogonal polynomials. In section 5 we compare 
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these approximate formulas with exact expressions calculated numerically for 
a = 1/3. We close the paper with a discussion. 



2 Universality regimes 

The universality of a correlation function in a random matrix model means 
that it does not change (up to a trivial rescaling) when the probability mea- 
sure is modified. Typically the universal quantities involve properties of 
the eigenvalues and their correlations considered on the scale of eigenvalue 
spacing in the large limit. For a multicritical point (see below), the modifi- 
cations of the probability measure are restricted to those that do not destroy 
multicriticality. However, in general, this still leaves an infinite dimensional 
space of possible deformations. 

The type of universal behaviour in random matrix models depends on 
the part of the spectrum that one is studying. In general the scaling of 
eigenvalue spacing with A^, in the vicinity of a fixed eigenvalue Aq, can be 
easily obtained by looking at the local behaviour of the (bulk) eigenvalue 
density (normalized to A^): 

N p{X) ^ N{X - XoT (1) 

Consequently the number of eigenvalues between Aq and A is approximately 
n ~ A^(A — Ao)""^^. Reexpressing A in terms of n shows that the eigenvalue 
spacing in the vicinity of Aq scales like l/N^^^^^^K 

In the bulk the eigenvalue spacing is of the order — this is the 
classical Wigner-Dyson regime ^ . At the edges of the spectrum the spacing 
is like l/N"^/^ and one observes universal Airy-like oscillations ^ [II]]- When 
the random matrix posseses an additional chiral structure as in applications 
to QCD i.e. 

there appears a different universal regime close to A = (with eigenvalue 
spacing ~ 1/A^) 0]. 

Apart from these universal regimes, when one finetunes the potential so 
that p(A) ~ A^™, a discrete series of multicritical points appear (with scaling 
\l]S[^/{'^'^+^)y The behaviour on the scale of eigenvalue spacing in this regime 
is very difficult to extract [H, IIU], |TT|. A different (but also discrete) class of 
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multicritical models appears when adding an appropriate fixed matrix to the 
random matrix. This behaviour and novel scaling has been analyzed both 
for ordinary |12[ and chiral 0] random matrix models. 



From the above discussion we see that in order to construct a class of 
random matrix models with arbitrary real scaling exponents we have to find 
a random matrix ensemble whose bulk eigenvalue spacing behaves like 

p(A) ~ |A|° (3) 

for real a. We will perform the construction for a G (0, 1) but an extension 
to other a should not be difficult. 



3 Eigenvalue density in the bulk 

For a hermitian random matrix model with measure 

^-NtTV(M) 

we consider the following class of potentials 

V{M) = --\M\^+^\M\^P (5) 
p 2p 

For hermitian matrices we may define \Ai\^ in the following way. By a 
unitary transformation Ai can be rewritten in the form Ai = U MJ'^ with A = 
diag{\i, . . . ,\n) a diagonal matrix. Then we may define \M.\^ = U\A\pW 
where |Ap = diag{\Xi\P, . . . , |A^rp). In fact we only need its trace which is 
given by 

N 

tT\M\P = J2\\,\P (6) 

i=l 

We note that although random matrix models with potentials of the type 
d^) are perfectly well defined, they seem to lack an evident diagrammatic 
perturbative expansion. Potentials with a single power-like term were first 
considered in the mathematical literature [0 and in The parameter p 
will be related to a by 



p=l+a (7) 

In the bulk we solve for the eigenvalue density by standard saddle point 
method. We are looking for a solution with a single cut {—a, a). Later we 
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will choose a = 1. The expression for the Green's function is then 

dX V'{X) 



G(z) = —Vz^ - a2 / , ^ ^ (8) 

^ ' 2n J-aZ-\ - A2 ^ ^ 

We note that for a > 0, the derivative exists. The cut endpoint is fixed 
through the constraint 

Trio - A2 ^ ^ 

The eigenvalue density can be extracted from the imaginary part through 

Because of the nonanalytic power like structure of the potential (^), it is 
convenient to get rid of the principal value and rewrite the formula as an 
ordinary integral. First using the fact that the potential is even we may use 



1 1 2A 



z-X z + X z^-X^ 

to obtain 



1 o r... r dX XV' (X) 

p(z) = -Va^ - z^ ■ PV — — , ^ ' 12) 



Now we use the fact that 



pyf ,.-0 (13) 

Jo z^ - A^ - A2 



Therefore we may rewrite (0) as 

. i^v^ r -'-'M - J(^) ^ / ^.^ (14) 

vr2 Jo z^ - A2 v^a2 - A2 

where we could erase the principal value as the integrand is now nonsingular. 

We will now fix the endpoint of the cut to a = 1. Substitution of (|) into 
d) yields 

— ^9 } .\ = 1 (15) 

5 



Since we want to locate a multicritical point we have to require that p(0) = 0. 
This gives the second equation for t and g: 

^'(^'''^ -0 (16) 
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I.e. 



The solution is 



Av^r^ 



rfe r tei 

t + g = (17) 

r(£) T{p) ^ ' 



(18) 
(19) 



2v/?pr(p-i) 


ir(f) 








1 - 2V 




)r(i±^) 



r(p + |)r(^)-2r(p-i)r(i±2) 



Bulk eigenvalue density 

Once the potential is fixed let us compute explicitly the bulk eigenvalue 
density 

p{z) = ^Vi - HM^) + gM^)] (20) 



TC 

where 



■py^) = /_ ... (21) 



/O Z2_p^/i^2 

This integral can be explicitly expressed in terms of hypergeometric functions 



M^) = )„\ 2i^i ( 1, 1 - ^ ) + o tan - — (22) 



2r 2 ^ ^ V' 2' 2 ' ; ' 2 2 ^-^ 



We see that it is the last term which gives the multicritical behaviour that 
we wanted to obtain p{z) ~ zP~^ = z". Owing to equation ([T6|) the constant 
terms cancel. We could have chosen of course a different power for the 
second term in (^, however with this choice the numerical construction of 
the relevant orthgonal polynomials is made easier. 

Let us note what happens when we reach the point a = 1. Then p = 2 
and the potential obtained here is the standard multicritical quartic one 
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V{x) = —Ax"^ + 4x^. However due to the vanishing of the coefficient of the 
last term in (^2]) the eigenvalue density behaves like p{z) ~ z"^ and not z. 



When we increase a the expression for the eigenvalue density ceases to be 
positive and thus the construction fails. Presumably a modification of the 
power in the second term of the potential might cure the problem but we 
will not consider that here. To sum up, the expression (pO]) is nonnegative 
when a G (0, 1) as considered in this paper. 

4 Behaviour near the origin — orthogonal 
polynomials 

The ultimate interest in constructing the random matrix models with the new 
scaling properties is to extract universal properties which typically occur in 
the microscopic regime i.e. on the scale of eigenvalue spacing. It is well known 
that all the relevant properties are encoded in the orthogonal polynomials 

/oo 
rfxP„(x)P™(x)e-^^(") = (23) 
-oo 



It is convenient to introduce the wavefunctions 

M^) = Pn(x)e-^^ (24) 



Then the kernel which allows for the determination of all correlation functions 
is given by the expression 

K[z,w) = \ Rn (25) 

* z — w 



where Rn is the recursion coefficient entering zipniz) = y/R^ipn+iiz) + 

^/R^i)n-l{z). 

The microscopic limit is defined through a rescaling 

X = -V (26) 

and a limit N oo with y kept fixed. The key quantity that determines the 
universal properties is the rescaled wavefunction 
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(as well as a similar expression with t/^tv-i), Cn is a normalization constant 
chosen in such a way that the limit exists. In fact the existence of such a 
limit is a nontrivial and nongeneric property. 

In general the determination of the scaling function F{y) is a very difficult 
problem c.f. ||, [l^, [ill. One can formulate a differential equation satisfied 
by the wavefunction tpni^)-, but taking the scaling limit is extremely difficult 
and e.g. for the case of quartic multicritical ensemble involves data from 
an auxiliary Painleve equation p, |10|. Moreover the starting point of such 
considerations requires a polynomial potential, which is not the case for our 
class of models. It would be interesting in the future to explore the possibility 



of applying the Riemann-Hilbert methods of [11 



Here we would like to adopt a different approach and give a very simple 
but approximate method of constructing the scaling function F{y). In section 
5 we will compare the approximate solution with the numerically obtained 
exact result for a = 1/3, and with the mesoscopic approximation of |^. 



Zeroes of the orthogonal polynomials 

The starting point of our construction is the elementary fact that we may 
reconstruct the orthogonal polynomial from the knowledge of its zeroes. Since 
in any case we normalize the (even) wave functions through il'2n{^) = 1 the 
appropiate formula isQ 



e 2 



n 1- 



1=1 



X 



(28) 



Interestingly enough there is a theorem due to Ismail [|15[ which states 
that for an (almost) arbitrary potential V{x) (see WW) the zeroes of the n*"^ 
orthogonal polynomial Aj solve the equations 



E 



NV'{\i 



l<k<n,i^k 



A,; — Al 



+ 2 ^^^^ 



(29) 



i.e. these are electrostatic equilibrium positions in an external field. The 
additional assumptions present in serve only to prove the uniqueness of 
a solution to (Bg|). 



^We always consider even potentials. 
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The correction term involves the function An{x), given by the formula 
AM^^r (30) 

"y -tlfi J —oo X — y 

For polynomial potentials of degree m, is a polynomial of at most 

degree m — 2. Note that the above equations are exact and valid for any A^. 
We see that the last term in (^) is suppressed w.r.t. the ordinary potential 
term. 

Now in order to develop our approximation scheme let us neglect the 
correction factor (logy4„(x))' (this is an approximation, the An term may 
indeed give some contributions in the scaling limit). From (^) we see that 
the zeroes are distributed with a continuum density identical to the eigen- 
value density p(A) of the corresponding random matrix model. The basis 
of our approximation scheme is the assumption that locally the zeroes are 
distributed uniformly w.r.t. p(A) (see the examples below). 

Classical Wigner-Dyson scaling functions 

Before we consider the case of immediate interest to us, in order to illustrate 
the method let us first rederive the formulas for the standard Wigner-Dyson 
scaling functions. To this end we assume that the eigenvalue density at A = 
is nonvanishing and is very large. For simplicity we take the potential to 
be even. It is convenient to normalize the continuum density of zeroes by 

Pcont 

{X)dX = N (31) 

In our approximation PcontW is simply given by ■ p(A) where p(A) is the 
bulk eigenvalue density of the random matrix model (normalized to 1). If 
the number of zeroes is even then they are distributed symmetrically around 
and there is no zero at A = 0. The assumption of 'uniform distribution' 
means that the i — th zero is located at A^ where 

Ai 1 

Pcont{X)dX = i + - (32) 

Since we are interested only in zeroes which lie close to A = 0, in the above 
expression we may substitute PcontW -^p(O)- Consequently Aj = (i -|- 
l/2)/(A^p(0)). Therefore the wavefunction is 



9 



Here we extended the upper limit of the product to infinity. After introducing 
the scahng variable y = xN we obtain immediately 



= cos(7rp(0)y) (34) 

the exact scaling function in the bulk. The odd case is similar but then there 
is a zero at A = and consequently the relevant equation is 

A,; 

p{X)d\ = i (35) 







From the infinite product representation one obtains then the sine function. 
The kernel (^5]) in the scaling limit follows immediately 



sin {7Tp{0)yi) cos {np{0)y2) - cos {np{0)yi) sin (7rp(0)?/2) 
K[yi,y2) ~ 



yi -y2 

sin(7rp(0)(yi - ^2)) 



(36) 



The fact that we obtained here the exact scaling functions does not mean, 
however, that we should expect to get exact results in general. 

Approximate scaling function for a E (0, 1) 

The analysis of the scaling function for the new multicritical ensembles is 
very similar. For definiteness we will just consider the even case. We start 
from the behaviour of the eigenvalue density near zero: 

p{X)=c^\Xr (37) 



where Cq, can be easily extracted from (pO]). The positions of the zeroes 



according to the 'uniform distribution' approximation are determined by 

/'^ Nc^X'^dX = -^c^NXf^' = -^c^yt^^ =i + ]- (38) 
JQ a + 1 a + 1 I 

where we introduced the rescaled variable y = XN^^^^'^°'\ Hence 

'1 + a / 1 M 



c„ V 2 
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Consequently the scaling function normalized by F{0) = 1 is 



F{y) = exp 



t y 



l+a ' 



1 + a 2 



1=0 \ yi J 



(40) 



In the next section we will compare this approximate result with the exact 
orthogonal polynomials for a = 1/3. However because the convergence prop- 
erties of the infinite product are not very good we will truncate it at some 
f^max (for numerical computation we use rimax = 200) and approximate the 
rest of the terms through 



The result is 



exp 



1 — a y„ ,1 



P(A)1 
(41) 

(42) 



As we see when taking the Umax oo limit this term becomes equal to unity 
(because yumax+i ~^ cxd). However we include it just because of the slow 
numerical convergence of the infinite product. 



5 An example 



a 



1 

3 



The orthogonal polynomials for any potential can be constructed using the 
determinant formula: 



(43) 





( Ao 


Ai 


A2 . 


■ A„ \ 






A2 


A3 . 


• A„+i 


~ det 


A2 


A3 


A4 . 


. A,,+2 








x"^ 





where A,- are the moments 



A, 



X e 



-NVix)^^ 



X 



(44) 



The advantage of the specific form of (^ is that the moments Aj can be 
expressed analytically in terms of the confluent hypergeometric functions iFi. 
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We constructed orthogonal polynomials for a = 1/3 (p = 4/3) for = 
40, 120 and = 240. In Fig. 1 we show the wavefunctions expressed in 
terms of the scaling variable y = xN^^^ 

rr^iy) = (y/iV^) = Pn {y/Ni) ■ e"^''('/^') (45) 

for the above values of A^, normalized hj iPn{0) = 1. We see convergence 
towards a well defined scaling function. 



-0 . 5 
-1 




Figure 1: The wavefunctions ''pffiy) for a = 1/3 expressed as a function of 
the scaling variable y for N = 40, 120 and = 240 (thick line). 

In Fig. 2 we compare the exact wavefunction for A^ = 240 with the 
approximate scaling function F{y) obtained in the previous section. There 
is no free parameter in F{y). The agreement is indeed surprisingly good. 
Even the small 'bump' close to ?/ = is correctly reproduced. In fact it is 
difficult to judge from the numerical comparision whether the approximation 
is exact or not in this case. The small deviations might be caused either by 
true corrections which go beyond our approximation or by finite size effects. 
The reason why such a simple approximation scheme works so well certainly 
deserves further study. 

Finally let us say a few words about the comparison with the 'mesoscopic 
approximation' of An analysis of the differential equations for ipNi^) for 
the potentials considered here is still lacking, so we will just take a suitable 
analytical continuation of the mesoscopic approximation from the discrete 
multicritical points considered in ^ to our case. The result is just 

Fmesoiy) = COS ( ^^C^y^+A (46) 

Va + 1 / 
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Figure 2: Comparison of the wavefunction '?/'24o(?/) (thick hne) with the ap- 
proximate scahng function F{y). 

We see that the zeroes of this function on the real axis coincide with the 
zeroes of our approximation, but the analytical structure in the complex 
plane is certainly different. There are spurious cuts and zeroes coming from 
?/^/^. Nevertheless this approximation is also quite good but the substructure 
close to ?/ = is not captured by the mesoscopic approximation (similarly 
as in the quartic multicritical case considered in |^). The comparison is 
presented in Fig. 3. 




Figure 3: Comparison of the wavefunction ip2%o{y) (thick line) with the meso- 
scopic scaling function 1.04 ■ Fmesoijj)- 
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6 Discussion 



In this paper we constructed a continuous class of random matrix ensembles 
which exhibit a new type of universal behaviour close to the zero eigenvalue. 
The eigenvalue spacing in this region behaves like l/A^^/("+^) with a real 
exponent a G (0,1). We derived the bulk eigenvalue density and studied 
approximately the behaviour of the scaling functions in the microscopic limit. 
Our method for obtaining the scaling function was approximate and based 
just on the local behaviour of the bulk eigenvalue distribution. Nevertheless 
the approximation seems to be very good and captures even the fine structure 
of the scaling functions near the origin. 

It would be very interesting to see if these new universality classes could 
appear in some physical systems. We conjecture that the scaling behaviour 
can be directly related to the large eigenvalue behaviour of Levy random 
matrix models [^. 

There are numerous open questions and directions for further study. 
Firstly it would be very interesting to try to obtain, even implicitly through 
e.g. a differential equation, the exact scaling function, especially as the vast 
majority of techniques dealing with multicritical random matrix ensembles 
requires from the outset a polynomial potential. Perhaps the most promis- 
ing approach would be the Riemann-Hilbert method. On the technical side 
a more systematic analysis of the properties of the orthogonal polynomial 
zeroes starting from (^), in particular the nature and scaling properties of 
corrections would be very welcome. This is especially interesting in order 
to better understand the effectiveness and limitiations of the approximate 
method which seems to work so well here. In general, however, this seems 
to be quite a formidable problem. Secondly other quantities of interest, such 
as eigenvalue spacing distributions, would be important in view of possible 
applications. Other possible directions are an extension of the above consid- 
erations to chiral random matrix models (in view of possible applications to 
QCD), and a rigorous treatment of universality properties. 

Note added: An independent study of a very similar class of multicrit- 
ical random matrix model with real exponents was undertaken by Gemot 



Akemann and Graziano Vernizzi |T^]. 
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